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Q . Abstract 

O 

^^ \ We present a detailed theoretical investigation of hadron attenuation in deep inelastic scatter- 

D ' ing (DIS) off complex nuclei in the kinematic regime of the HERMES experiment. The analysis 

pq ■ is carried out in the framework of a probabilistic coupled-channel transport model based on the 

Boltzmann-Uehling-Uhlenbeck (BUU) equation, which allows for a treatment of the final-state 

(N 

K^ ' interactions (FSI) beyond simple absorption mechanisms. Furthermore, our event-by-event simu- 

m ■ 

^^ ' lations account for the kinematic cuts of the experiments as well as the geometrical acceptance of 

O 

^t \ relative to deuterium as a function of the photon energy u, the hadron energy fraction z^ = Ey^jv 



the detectors. We calculate the multiplicity ratios of charged hadrons for various nuclear targets 



and the transverse momentum p^. We also confront our model results on double-hadron attenua- 

^ ■ tion with recent experimental data. Separately, we compare the attenuation of identified hadrons 

5^ . (tt^, vr*^, i^^, "p and "p) on ^'^Ne and ^^Kr targets with the data from the HERMES Collaboration 

^_, and make predictions for a ^^^Xe target. At the end we turn towards hadron attenuation on ^^Cu 

^ • nuclei at EMC energies. Our studies demonstrate that (pre-)hadronic final-state interactions play 

C^ " 

a dominant role in the kinematic regime of the HERMES experiment while our present approach 

overestimates the attenuation at EMC energies. 
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I. INTRODUCTION 

During the past years the HERMES collaboration has carried out a detailed experimental 
investigation of deep inelastic lepton scattering off complex nuclei [l|, ^ . This has lead to 
numerous excellent data on the attenuation of high-energy hadrons in 'cold' nuclear matter 
which can be used to study the hadronization of particles created in (hard) collisions as 
proposed in Ref. J3j. 

Besides being an interesting topic on its own the detailed understanding of the space- 
time picture of hadronization should also help to clarify to what extent the jet suppression 
observed in heavy-ion collisions at RHIC M is due to (pre-)hadronic final-state interactions 
(FSI) p, |6| or partonic energy loss [Tj . In this respect the lepton scattering experiments also 
provide a benchmark for the interpretation of ultra-relativistic nucleus-nucleus collisions 
that aim at probing the transition from deconfined matter to a color neutral interacting 
hadron gas. 

Contrary to nucleus-nucleus collisions, where the created 'fireball' is rapidly expanding, 
lepton-nucleus reactions provide rather clear geometrical constraints, which are well un- 
der control experimentally. In the one-photon exchange approximation the scattering of a 
charged lepton off a nucleon or nucleus can be decomposed into the emission and absorption 
of a virtual photon with energy z/, virtuality Q^ and polarization e. In such a high-energy 
photon-nucleon interaction the reaction products hadronize, due to confinement, long before 
they reach the detector. One can estimate the formation time of a hadron by the time that 
its partonic constituents need to travel a distance of the order of a hadronic radius. Ac- 
cording to this estimate the formation time in the rest frame of the hadron is about 0.5-0.8 
fm/c. Because of time dilatation, the formation length in the laboratory frame then can ex- 
ceed nuclear radii at high energies. In lepton-nucleus interactions the reaction products can, 
therefore, interact with the surrounding nuclear medium during the formation time. Here 
the nuclear target can be viewed as a micro-laboratory that provides an intrinsic (variable) 
time scale due to the size of the target nuclei, which can be exploited to get information on 
the actual time scale of the hadronization itself. 

The HERMES collaboration has measured the multiplicity ratios of pions, charged kaons, 
protons and antiprotons in deep inelastic positron scattering at 27.6 GeV off ^D, ^^N, ^°Ne 
and ^^Kr targets. The observed attenuation of the hadron multiplicities in complex nuclei 



has lead to different interpretations. One possible explanation is that the quark, that was 
struck by the photon, propagates through the nucleus and undergoes multiple scattering 
[8|, |9[, thereby loosing energy by induced gluon radiation. The authors of Ref. [10| attribute 

hadron absorption and a rescaling of the 
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we have shown that the attenuation 



the observed attenuation to a combined effect o 

quark fragmentation function in nuclei. In Refs. 

of charged hadrons at HERMES energies can also be understood by final-state interactions 

(FSI) of color neutral prehadrons that are produced quite early after the photon-nucleon 

interaction. This very short production time is in line with the gluon-bremsstrahlung model 

of Ref. [13. 

In this work we discuss different concepts of hadron production and formation in the 
framework of a probabilistic coupled-channel transport model which allows for a realistic 
treatment of the FSI. We demonstrate that the coupled-channel effects in the FSI lead 
to deviations from the results of simple absorption models since particles produced in the 
FSI affect the low energy part of the hadron spectra and thereby all energy-integrated 
observables. A strong effect of the coupled channels becomes also visible in the attenuation 
of kaons since kaon absorption is partially compensated by kaon recreation in the FSI via 
secondary channels like tc + N —^ K + Y + X. Furthermore, our event-by-event simulation 
allows us to investigate how the observed attenuation ratios depend on the kinematic cuts 
and the detector acceptance. We will show that there is a sizable effect of the detector 
acceptance on the multiplicity ratios, which does not drop out by taking the multiplicity ratio 
and, therefore, has to be considered explicitly when addressing any robust interpretation of 
the data. 

Being able to compare with almost every observable that is experimentally accessible we 
can test the limits of a purely hadronic model and figure out at what point an extension 



becomes necessary. For this reason we also include the EMC data [14] on hadron attenuation 
at 100 and 200 GeV muon energy in our analysis where the hard scale - set by the photon 
virtuality Q^ - is closer to that of the high pt events at RHIC 4]. We stress that our model 
is also applicable for a broad range of other reactions such as vrA, pA and AA collisions. A 
related study of hadron attenuation for high transverse momentum particles m. d + Au and 
Au + Au collisions at RHIC energies has been carried out in J6[ . 

Our paper is structured in the following way. In Sec. |H] we discuss the description of 
the initial (virtual) photon-nucleon interaction using the event generators PYTHIA 6.2 and 



FRITIOF 7.02 and test our model input with experimental HERMES data on a hydrogen 
target [15[ . Section IIIII includes a brief description of the shadowing corrections employed 
for reactions on nuclei. The BUU transport model - used for the description of coupled 
channel FSI - is recalled in Sec. IIVI while in Sec. |V] we discuss our concept of prehadrons. 
In Sec. IVll and Sec. IVllI we compare the results of our model with experimental data of 
the HERMES and EMC collaborations. We focus on the multiplicity spectra of charged 
hadrons as a function of the photon energy z/, the hadron energy fraction Zh = Eh/v and the 
transverse momentum px- Furthermore, we discuss the effect of the detector geometry and 
investigate the double-hadron attenuation as well as the attenuation of identified hadrons 
such as TT^, vr*^, K^^ p and p. We close with a short summary in Sec. IVlllI More technical 
details of the transport approach are presented in the appendices. 

II. HIGH-ENERGY PHOTON-NUCLEON INTERACTIONS 



Based on the ear 
oped in Refs. 
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ier work on high-energy photon-nucleus reactions in [Xol we have devel- 
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19| a method to describe high-energy photon- and electron-induced 
reactions in the framework of a semi-classical transport model. We thereby split the inter- 
action of the virtual photon with the nucleus into two parts. In the first step i) the virtual 
photon interacts with a bound nucleon of the target to produce a final state which in step 
ii) is propagated within the transport model. In the initial 7*A^ interaction we account for 
nuclear effects, such as Fermi motion, Pauli blocking and binding energies. Furthermore, 
coherence length effects in the entrance channel - leading to nuclear shadowing - are taken 



care of using the method developed in Refs. 



la 



m- 



For small invariant masses of the photon- nucleon system {W < 2 GeV) we use an explicit 
resonance description for the primary (virtual) photon-nucleon interaction as in Ref. J20|. 
Above the resonance region 7*A^ interactions become more complex and in our model the 
final state is in general determined by the event generator PYTHIA 6.2 21]. We recall that 
in a high-energy 7*A^ interaction the virtual photon does not always couple directly to a 
quark of the target nucleon; depending on its kinematics it may also fluctuate into a vector 
meson V = p^,u!,(f), J/tp (vector meson dominance (VMD)) or perturbatively branch into 
a quark-antiquark pair (generalized vector meson dominance (GVMD)) which subsequently 
scatters off the nucleon. In case of VMD events the same processes as in hadron-hadron 



interactions can occur including diffractive scattering and hard interactions between the con- 
stituents of the nucleon and the vector meson. In addition, most of the hard interactions in 
PYTHIA involve gluon-bremsstrahlung processes. The high-energy interactions of the pho- 
ton, therefore, lead to the excitation of several hadronic strings with large invariant masses 



22J | routine JETSET implemented 



whose decays are described by the Lund fragmentation 
in PYTHIA. 

For technical reasons, resolved 7*A^ events, i.e. VMD and GVMD events, can only be 
generated by PYTHIA if the invariant mass W of the photon-nucleon system is larger than 
a threshold energy Wpy which depends on the PYTHIA parameters; it takes the value 
WpY = 4 GeV for the default parameter set. For invariant masses 2 GeV < W < Wpy we 
employ the event generator FRITIOF 7.02 [23] to simulate the resolved photon interactions. 
As discussed in Appendix |X] FRITIOF is also used to model the high-energy hadronic FSI 
in case of nuclear reactions. In contrast to PYTHIA, the FRITIOF model can be applied 
down to the onset of the resonance region. Below invariant energies of 10 GeV the hadron- 
hadron interactions in FRITIOF are alwa ys d iffractive. Throughout this work we employ 



the same FRITIOF parameters as in Refs. yA 



2^. InRef. M 



In Ref. 241 it was shown that this choice 



of parameters provides a good description of particle production over a broad energy range. 
The only exception is the suppression factor for ss creation in the string fragmentation, 
which we reset to the FRITIOF default value 0.3 to improve the agreement with more 
recent strangeness production data |26l |. 

To simulate the VMD part of the 7*A^ interaction with FRITIOF we pass the photon as 
a vector meson V with a probability 

where the numerator, according to Refs. J2l|, |27|, gives the contribution of the vector meson 



component V to the total photon-nucleon cross section 0"*°*jy. In Eq. ([T)) the VMD coup 



constants are denoted as gy] o"y^ are the total vector meson nucleon cross sections 
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mg 



^poV(s) ~ al°^is) ^ 13.63/ + 31.79s-'' [mb], 
(tI%{s) ^ 10.01s' - 1.52s"'' [mb], 

^"/W(^) ^ 1^C(^)' (2) 



with s in GeV^ and e = 0.0808, r] = 0.4525. The factor 

{ml. + Q2)2 



rK,Q^)=0.5 ,J, 1^^,,, (3) 



accounts for the contribution of longitudinal photons to the photon-nucleon cross section. 

In Appendix^ we show that FRITIOF insufficiently describes elastic scattering, which 
directly affects the diffractive vector meson photoproduction 7A^ —>■ VN. In the transition 
region below Wpy "we, therefore, parameterize the diffractive production cross section for p^, 
ijj and (\) using the Regge prescription of Donnachie and Landshoff 28| which also provides the 
correct dependence on the four-momentum transfer squared t. The high-energy cross sections 
are continuously extrapolated to our parameterization in the resonance region (cf. Ref. (iq). 
For the description of J/-?/' production close to threshold we use the two and three gluon 



exchange model of Brodsky et al. 



and parameterize the high-energy y^ dependence as 



0.77 



cT^7v^j/^7v = 0.002 (^^j /ib. (4) 

Fig. n shows our parameterization of diffractive vector meson production for p°, tu, and 
J/?/' in comparison to a collection of experimental data and demonstrates that our 'input' 
for the 7*p reaction is well in accordance with experimental information. 

For the processes 7A^ -^ V IS. we use our parameterizations of the differential cross section 
for 7A^ — > 0A and 7A^ -^ cjA from Ref. [3l| and take the same matrix element for 7A^ -^ p" A 
as for 7A^ — » tuA. 

To account for the virtuality Q^ of the photon we multiply the vector meson photoproduc- 
tion cross sections with the VMD formfactor of Ref. J27|, which is also used in the PYTHIA 
model (cf. Eq. (H))) 

The FRITIOF model has the tendency to simply flip the spin of the incoming vector 
meson to zero thereby creating a lot of events like 7*A^ -^ tt^N with vanishing momentum 
transfer, a problem not present in PYTHIA. These events would correspond to diffractive 
7r° production. However, since the Pomeron carries the quantum numbers of the vacuum 
it cannot change the charge conjugation quantum number of the incoming p'' (C = — 1) to 
that of the vr'^ (C = +1). We, therefore, simply remove these unphysical events from our 
simulation. 



We have compared the particle spectra of VMD events as generated by FRITIOF and 
PYTHIA at the threshold energy Wpy = 4 GeV and found no significant differences between 
the two approaches. However, we keep the option to also generate the VMD events above 
WpY = 4 GeV by FRITIOF to ensure a continuous transition at invariant energies W > Wpy 
For calculations at EMC energies we always use the PYTHIA model for event generation. 
Note, however, that the VMD contribution vanishes with increasing Q^ and already at 
Q^ = 1 GeV^ the VMD part aZ-^^ -n accounts for only 50% of the total photon-nucleon cross 



section at HERMES energies 19 1. 



The GVMD part below Wpy is modeled by passing a vector meson with the quark content 
of the qq fluctuation to FRITIOF. In analogy to Eq. ^ the probability is now given by the 
corresponding GVMD term of Ref. 27[ normalized to cr*°*^. As discussed in Ref. |22l] the 
products of a GVMD and VMD event are different in nature since the VMD and GVMD 
components of the photon involve different intrinsic transverse momenta. However, we 
have shown in Ref. [l3] that the GVMD part in the transition region only contributes by 
less than 15% to the total photon-nucleon cross section. Hence, one can postpone a more 
sophisticated treatment unless one concentrates on events that are especially triggered by 
the GVMD component of the photon. 

We now compare our results for electroproduction on a proton target with the hydrogen 
data of the HERMES collaboration [l5| taken at a positron beam energy -Ebeam = 27.6 GeV. 
In Fig. 121 we show the Zh = Eh/v spectra dNh/dzh of n^, K^, p and p normalized to the 
number Ne of deep inelastically scattered positrons. In our calculation we apply the same 
cuts on the event kinematics as in the experiment: For the positrons these are Q^ > 1 
GeV^, W'^ > 10 GeV^ and 0.1 < y = ///-Ebeam < 0.85. For the hadrons we require 1 GeV/c 
< Pt, < 15 GeV/c, 2 GeV/c < PK,p,p < 15 GeV/c as well as xp > 0.1, where we define the 
Feynman variable xp as in experiment by 

pT 

xf = -p — I ■ (6) 

IQ'cm I 

Here py™ denotes the momentum of the hadron parallel to the momentum g^m of the virtual 
photon in the center-of-mass frame of the photon-nucleon system. Since the data are not 
acceptance and efficiency corrected we account for the angular acceptance of the HERMES 
detector j32.], i.e. ±170 mrad horizontally and ±(40 — 140) mrad vertically, for both the 
scattered positrons and the produced hadrons in our simulation. The solid line in Fig. |21 



shows our result using FRITIOF for the resolved events below Wpy = 4 GeV as well as for 
all VMD events above Wpy- The dashed line represents a calculation where we have changed 
the default parameters of PYTHIA 6.2 in such a way that it is applicable down to Wpy ~ 3 
GeV and, hence, simulate all events with PYTHIA. Except for a small deviation in the 
proton spectra both methods yield essentially the same result. In view of the fact, that we 
do not include the detector efficiency in our calculation, which is unknown to the authors, but 
only account for its angular acceptance, our calculations are in satisfying agreement with 
the experimental data. One can see that without any further fine tuning our simulation 
including the kinematic cuts and detector acceptance reproduces the absolute size of the 
multiplicity spectra. We have also compared the xp, pt and total momentum-spectra of 
the different hadron species with the HERMES data and find a similar good agreement. 
To demonstrate the effects of the kinematic cuts and the limited angular acceptance of the 
HERMES detector, the dotted line in Fig. |21 shows the results of a simulation where no cuts 
on the hadron kinematics have been applied and where we have assumed a 47r-detector. 

III. NUCLEAR SHADOWING 

As discussed in the previous section the photon does not always interact directly with a 
quark inside the target nucleon since its wavefunction is a superposition of a bare photon 
and resolved hadronic fluctuations: 

I7) =Cbare|7bare)+ JZ Cy|\/)+ ^ Cg\qq) . (7) 

V=p'^,ijj,<j),J/'tl' g=u,d,s,c,b 

In an electron-nucleus interaction the hadronic components of the photon in ((7j) are shadowed 
on their way through the nucleus due to the strong interaction of these components with 
the hadronic environment. The strength of this shadowing effect depends on the 'coherence 
length', i.e. the distance Ih that the photon travels as a hadronic fluctuation h. This distance 
can be estimated via the uncertainty principle: 

Ih = Qh' = {k - h)-\ (8) 

where k = ^v"^ + Q^ is the photon momentum and k^ = \/v'^ — "m^ denotes the momentum 
of the hadronic fluctuation on its mass shell. If the coherence length l^ becomes larger 
than the mean-free-path of the hadronic fluctuation in the nuclear medium the hadronic 
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fluctuation gets shadowed. In the kinematic region of interest in this work the shadowing 
of the massive qq fluctuations can be ne glec ted and one is only left with the shadowing of 
the photon's VMD components. In Ref. jlSl . |l9| we have calculated the modification of the 
vector meson components at position f = {b, z) in the nucleus: 



\V)^(l-T^^\f^)\V), 
where the nuclear profile function Ty is determined within Glauber theory, i.e. 



^(A) 



T'y\h,z) 



dzip{b, Zi 



a 



tot 

VN , 



X exp 



--<J 



_tot 



iav)e 



iqv(zi-z) 



VN 



(1 - iay) 



z 

/ dzkp{h, 



Zk) 



(9) 



Here a^^ denotes the total vector meson-nucleon cross section and ay is the ratio of the 
real and imaginary part of the VN forward scattering amplitude. Note that the coherence 
length (jSJ enters in the phase factor of Eq. 0. The modification of the photon's vector 
meson components are then taken into account when genera ting the scattering event with 
PYTHIA and FRITIOF. As we have demonstrated in Ref. [19| our prescription is in full 
agreement with the coherence length effects observed in incoherent p^ production off nuclei. 



IV. BUU TRANSPORT MODEL 



Our transport model - employed for the description of FSI - is based on the Boltzmann- 
Uehling-Uhlenbeck (BUU) equation. It has been originally developed to model heavy-ion 



Q 



collisions at SIS energies j^ and later been extended to describe pion-induced reactions 
as well as photo- and electroproduction J2i| in the resonance region. In the more recent past 
we have extended the model to simulate also nuclear interactions of high-energy photons and 
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Q 



electrons in the kinematic regime of the Jefferson Lab [16|, [Tfl, |l8| and HERMES 
experiments. Obviously, one of the advantages of the model is its applicability to many 
different nuclear reactions with the same set of parameters and physics assumptions. In 



addition, the coupled-channel treatment of the FSI goes far beyond the standard Glauber 
approach since it allows for a side-feeding of channels under study and not just absorption. 



The BUU transport model is based on a set of generalized transport equations for each 
particle species i, 

r^ + Vj^HVr- - Vr-HVA Fi{f, p, /i; t) = Icoll{{Fj}) (10) 

where 

H=^{fi + Us{f,p;t)y+fP (11) 

denotes the relativistic Hamilton function of a particle with mass /i in a scalar potential 
Us- For vanishing collision term Icon Eq. (jiup describes the time evolution of the spectral 
phase-space density 

Fi (r , p,fi;t) = fi (f , p; t)Ai (/i, p) (12) 

of non- interacting particles that move in a scalar mean- field potential Us, which in general 
depends on the phase-space densities {fj} of all other particle species (including i). The 
spectral function Ai of particle type i is parameterized as in Ref. |l6l |: 

A ( N _ 2 /i^rtot(/i,p) . . 

where Mj denotes the pole mass of particle i and Ftot its total width which in medium also 
depends on the particle momentum p. 
The collision term 

/coii({i^,}) = -t^fi{F,})Ah - 2S<({F,})A(1 ± h) (14) 

in Eq. (fTUj) accounts for changes in the spectral phase-space density of particle species i = ai 
due to multi-particle collisions of the type ai, . . . , a^ — > 6i . . . 6m and bi . . .km ^ ai, . . . ,an. 
The factor (1 ± /) has the plus sign for bosons (Bose enhancement) and the minus sign for 
fermions (Pauli blocking). The collision term consists of a loss term for particle species ai 

n m 

x(2vr)V(^Pa^. - ^P6j^a2,...,a„^bi,...,fe™|-Mai,...,a„^6i...fe™P 
i=l fc=l 
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and a gain term 




- ^K = 


1 


2Ea, 



f] f) cj i I 

\k=2 ^ ' "-k 



3=1 k=l 



12 
+ai...a„ I 



where gj are degeneration factors that account for the spin of particle j. The quantity |7Wp 
is the squared matrix element averaged over incoming and summed over outgoing spins. S 
denotes symmetry factors that take into account the number of identical particles in the 
incoming and outgoing channel, e.g., 

m ^ 

Sb^,...,b^ = n jri (17) 

k=l "*■ 
with Mb^ denoting the multiplicity of particle 6^. 

All transport equations (jTUj) are coupled by the collision term and the scalar potential 
Us in the Hamilton function (jllll where in our model the latter is only incorporated in case 
of baryons. Furthermore, in the general collision term we will restrict to 2 — > n transitions 
{n > 1) which is sufficient for the lepton- nucleus reactions of interest. For a specification 
of the scalar potential Us in case of baryons as well the 2 —^ n transition rates we refer the 
reader to Appendix El of this work. 



rehadrons' in the Lund model 




2^ with 



V. PREHADRONS 

In this section we compare the concept o: 
that generally used in transport models |l7l . l25l . |35[. The two concepts mainly differ in 
the treatment of the production proper time Tp of the color neutral prehadrons which are 
created in the string decay. While the production time of 'prehadrons' in the Lund model 
depends on the energy and momentum of the fragments, Tp is set to zero in conventional 
transport approaches |l7|, |25|, |35[- In order to specify the differences and to define the 



relevant quantities (production and formation times) we give a brief reminder of the Lund 
model in the following. 
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A. Prehadrons in the Lund model 



According to the Lund model [23] the hadronic strings - created in the high-energy 
interactions - decay due to the creation of quark-antiquark pairs from the vacuum. Fig. El 
shows the space-time picture of a Lund fragmentation process in the rest frame of a goQ'o 
string. In our example we neglect the masses of the quarks, i.e. all quarks move on the 
lightcone. At time zero the quark go carries the lightcone momentum p^ while the antiquark 
go carries the lightcone momentum p^. During their separation they loose energy and 
momentum to the string until they reach their turning points: 



4 = — , ^o' = (18) 



and 



< = ^. ^„- = ^, (19) 

where the string tension k in vacuum is known to be k ~ 1 GeV/fm. 

The production process of the n mesons with masses mi, . . . ,mn due to the creation of 
n-1 quark-antiquark pairs QjQj at the production vertices x • can be viewed as a series of 
steps along the (positive) hghtcone. Starting at the turning point (fT^ of the original quark 
go one takes random steps along the positive lightcone 

Axl = z,x+_i , z, e [0, 1] (20) 

where the probability distribution for the random variable Zj is given by the Lund fragmen- 
tation function 

/(^,)=Ari^^^exp(-6^^), (21) 

which depends on the transverse mass 

rriTj = Jrnj + p\ (22) 

of the produced hadron. A step along the negative lightcone is fixed by the on-shell condition 

for hadron j 

2 

777™. 

(23) 



(24) 
(25) 



AxJ = 


ml- 

K^Axj ■ 


Thus we obtain the recursion formulae 




4 = ^u - ^^/ 


= (1 - ^M-l 


X] = xj_i - AxJ 


-X- 1 ^-^^•'"^^■ 



12 



together with the constraint that the path must end at the turning point ()19|) of go- 
Obviously, the proper times rp{j) of the QjQj production vertices 

T^U) = x^xj (26) 

strongly depend on the mass rrij of the produced hadron. The same holds true for the proper 
time of crossing world lines for the pair qj_i and qj 

Tf{j)=xl_,xj. (27) 

The latter is often identified with the formation time Tf{j) of the hadron j. 

In the Lund picture there are in principle three timescales involved in the production of 
the hadron {qjqj_i). This is the formation time Tf{j) of the hadron as well as the times 
Tp{j — 1) and Tp{j) when the first and the second constituent parton is produced. As soon 
as the qj-i and the qj are created they form an intermediate object which can be viewed 



36, 



the longitudinal dimension of 



as a color neutral prehadron. As pointed out in Refs. 
the string and its (color neutral) fragments is never larger than typically 1 fm. Therefore, it 
does not make much sense to talk about freely propagating colored quarks between Tp{j — 1) 
and Tp{j) since each color charge is always attached via a string to the nearby anticharge. 
Instead one has to deal with the propagation of color neutral prehadrons and color neutral 
remainder strings during the fragmentation process. 

As an example we recall the production and formation times of the first rank particle 
(containing the go)) 

(28) 

(29) 

where for the first rank particle z corresponds to the fraction 2;^ of the positive lightcone 
momentum Pq that is carried away by the first rank hadron. For fixed z both times increase 
linearly with the (transverse) hadron mass rriT which implies that heavier hadrons are in 
general created later in time. Note, however, that the mass also enters into the Lund 
fragmentation function (j2^- Furthermore, the average production and formation proper 
times increase with the rank of the hadron as long as one considers a fragmentation into 
infinitely many hadrons. We stress that the assumption of infinitely many fragments enters 
explicitly in the derivation of the production and formation times in Ref. 38[ which are 
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r^ - 


m^ 1 — z 


'p 


K^ z 

m% 1 

Kj z 



averaged over infinitely many ranks. At HERMES energies, however, a string typically 
fragments into only 3-5 hadrons. Thus, when using the Lund production and formation 
times in our simulation, we explicitly extract the space-time points of the vertices and 
crossing points of pairs from the JETSET routines for each fragmenting string since the 
limit addressed in Ref. j381] is inappropriate for the kinematical regime of interest. 

In the upper panel of Fig. El we show the Zh dependence of the extracted average pro- 
duction proper times of the two hadron constituents as well as the hadron formation proper 
time Tf in electron-nucleon scattering at HERMES energies. We always label the smaller 
production time Tpi and the large Tp2. Note that, if more than one string is produced in 
the electron-nucleon interaction, the observed Zh = Eh/u is in general different from 2;^. 
The difference is most significant for z^ ^ 0.5. However, one also has to be careful with 
the interpretation of Zh in the limit Zh —^ 1. In this kinematic regime the hadron spectrum 
is dominated by diffractively produced vector mesons (cf. Fig. El) whose production times 
Tpi and Tp2 for the constituents are both zero. Furthermore, we have learned from our the- 
oretical study of incoherent p" electroproduction [19] that also the formation time Tf of a 
diffractively produced vector meson is compatible with zero in the HERMES kinematical 
regime. Consequently, the proper times Tpi, Tp2 and Tf for the p° (dash-dot-dotted line) 
vanish as ZpO ^ 1. 

In a direct photon-nucleon interaction the probability, that a meson contains the struck 
quark from the string end, increases with z^ and, consequently, Tpi vanishes for larger z^. 
On the other hand, the diquark that is left behind, i.e. the second string end, will most likely 
form a proton with small Zh which leads to a vanishing of Tpi for low-Zh protons (dash-dotted 
line). In case of diffractive photon-nucleon interactions the proton stays intact and Tpi, Tp2 
and Tf are zero which leads to a decrease of Tp2 and Tf for low-z/i protons. The Lund model 
predicts a vanishing of the production time Tp2 for z —* 1 which can be seen from Eq. ()28p. 
However, because of the finite bin-size the production time Tp2 in Fig. HI stays finite in the 
last Zh bin for all hadrons except for the diffractively produced vector mesons. 

In fact, one observes the general tendency that more massive hadrons like protons and 
p mesons are in general formed later in proper time than lighter hadrons like pions (solid 
line). This was already indicated by Eqs. ()281) and (j^^ for the production and formation 
times of the first rank hadron in a string fragmentation. Again the diffractively produced 
vector mesons are the exception. 
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We find that the production proper time Tpi of the first hadron constituent - averaged 
over all hadrons - lies between 0.04 fm/c at Zh ~ 0.95 and 0.4 fm/c at Zh ~ 0.35. The 
production proper time Tp2 of the second constituent is somewhat larger. It ranges from 0.3 
fm/c at Zh ~ 0.95 to about 1.2 fm/c at z^ ~ 0.35. The average hadron formation proper time 
Tf is of the order 1.1-1.5 fm/c except for the very large-Zh hadrons which are dominated by 
diffractively produced p° mesons and, hence, have again a very small formation time around 
0.4 fm/c at Zh ^ 0.95. 

In the lower panel of Fig. ^ we also show the corresponding production and formation 
times in the laboratory frame, i.e. the rest frame of the target nucleon: 

ip = 1-Tp = Tp (30) 

nih 

tf = -i ■Tf = Tf. (31) 

' ' mh 

We emphasize that simply boosting the proper times with the Lorentz gamma of the hadron 
- as done in Eqs. (IHUj) and (jHH) - is just an approximation. We currently carry out a 
more thorough investigation on this topic J39|. Because of time dilatation the light pions 
now have average production and formation times as large as 10-70 fm/c. However, as we 
discuss in the next section (and as can be seen from Fig. ^ pion production receives a large 
contribution from string fragmentation into p mesons that - due to their larger mass - have 
a considerably smaller average production time 0.5-3 fm/c in the laboratory frame. 

Note that the production and formation times in Fig. |3] are averaged iimes. In our actual 
simulation we assign in each scattering event the individual production and formation time 
from JETSET for each hadron. 

B. Prehadrons in transport models 

The hadrons containing the endpoint partons go and go of the string, that emerges from 
the excitation of a hadron in a scattering event, are not entirely created from the vacuum 
as the other qq pairs. Indeed, the hadrons that emerge from the string ends can be viewed 
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35| each string 



as the remnants of the original hadron. In the transport models |lj 
decay into color neutral prehadrons is determined instantaneously at the production vertex. 
During the formation time Tf only the beam and target remnants are allowed to interact 
with the nuclear environment whereas the hadrons containing only quarks or antiquarks 
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from the vacuum are propagated freely without interactions up to their hadron formation 
time in the computational reference frame. 

It is hard to identify the beam and target remnants of a 7*A^ event because particle 
production in PYTHIA and FRITIOF is complicated by the resolved photon interactions as 
well as the initial and final state gluon radiation that occur in many hard scattering events. 
As a result one often ends up with more than a single string per DIS event or a string that 
contains one or several radiated gluons, i.e. it might have the structure qgq, qggq etc. In 
general, the gluons have momentum components transverse to the cm momentum of the 
quark and antiquark. This leads to more complicated string topologies than just a linearly 



expanding string [2^ and has to be taken into account in the fragmentation. 

To identify the target and beam remnants of a binary collision we, therefore, trace the 
(anti-)quarks of the projectile and target all the way through the fragmentation process in 
JETSET. At the end of the reaction we then know those hadrons that contain the original 
(anti-)quarks, i.e. the beam and target remnants. Correspondingly, also those hadrons 
exclusively made from (anti-)quarks created in the string fragmentation are known explicitly. 

In Fig.Elwe show the energy spectra of hadrons produced in electron-nucleon interactions 
at HERMES. The different lines indicate the hadrons that contain zero (solid line), one 
(dashed line), two (dotted line) or three quarks (dash-dotted line) from the beam or target. 
In the dotted curves one can clearly identify the diffractive peak in the p" and (p spectra 
at z/t ~ 1, where the vector mesons contain the qq of the resolved photon. By comparing 
the p° and 7r° spectra one finds that the p° strongly contributes to the total pion yield by 
its subsequent decay into tt^tt^. From the K~^ and K~ spectra one sees (dashed lines) 
that due to their quark content (us) the K~ mesons contain less quarks from the beam or 
target than the K^ mesons. The K~, that are not solely made of quarks and antiquarks 
created from the vacuum in the string fragmentation, carry (anti-) quarks from the resolved 
photon component or the nucleon sea. Finally, one finds that there are only very few protons 
containing no quarks from the beam or target since diquark-antidiquark creation is strongly 
suppressed in the string fragmentation due to the relatively large diquark masses. This also 
explains why most protons at large Zh contain two of the original quarks, i.e. the diquark 
from the target nucleon struck by the photon. 

In our default approach we set the production times Tp of all prehadrons to zero and 
rescale their cross sections during the formation time r/ according to the constituent quark 
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model: 



tot _ ^org 

^prebaryon ~ o '-'baryon. 



3 

^orj 

T 



tot _ ' "org (on\ 

'-'premeson ~ r, '-' meson ) {'^■^J 



where norg denotes the number of (anti-) quarks in the prehadron stemming from the beam or 
target. As a consequence the prehadrons that solely contain (anti-)quarks produced from the 
vacuum in the string fragmentation do not interact during Tf. The assumption (jH^ is guided 
by the constraints of unitarity, which implies that the summed cross section of the products 
from a scattering process should not exceed the cross section for the initial hadrons on short 
time scales. For simplicity we assume that the formation time is a constant Tf in the rest 
frame of each hadron and that it does not depend on the particle species as in Refs. 2J, |25| . 



In Section IVllI we also discuss other concepts of the prehadronic cross section. However, 
as we will show in Sec. IVII the 'standard' concept (J32|l already gives a good description of 
the available HERMES data. This also holds for particle production in proton-nucleus and 
nucleus- nucleus collisions in a wide dynamical range j4fl( |. 

VI. HADRON ATTENUATION AT HERMES ENERGIES 

We start with an investigation of charged hadron attenuation in semi-inclusive DIS of 
27.6 GeV positrons off nitrogen and krypton. We apply the same kinematic cuts as in 
experiment 1^: Q^ > 1 GeV^ W > 2 GeV, u > 7 GeV, y < 0.85, x > 0.02 and Eh > 1.4 
GeV as well as the cut Zh > 0.2 for the u and pt spectra. Furthermore, we account for the 
angular acceptance of the HERMES detector. The left panel in Fig. IHl shows the average 
values of the photon energy u and the virtuality Q^ as a function of the energy fraction Zh 
of the charged hadrons produced on a ^^Kr target; the right panel shows the average values 
of Zh and Q^ as a function of u. In the simulation we used the prehadron concept (jH^ and 
a formation time Tf = 0.5 fm/c. Obviously, our simulation is in perfect agreement with the 
experimental values for the average kinematic variables (open symbols). 

The observable of interest is the multiplicity ratio defined as: 

Tjh / 2 /^2\ Ne(u,Q'^) \a /oo\ 

R^,{zh,iy,PT,Q ) = ^^(,^,.,,^,g.) . (33) 

Ne{u,Q^) Id 
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where A^^ is the yield of semi-inclusive hadrons in a given {zh, ViPr-, Q^)-bin and A^e the yield 
of inclusive deep inelastic scattering leptons in the same (z/, Q^)-bin. For the deuterium 
target, i.e. the nominator of Eq. (jHH|) . we simply use the isospin averaged results of a proton 
and a neutron target. Thus in the case of deuterium we neglect the FSl of the produced 
hadrons and also the effect of shadowing and Fermi motion. 

A. No prehadronic interactions 

Before discussing the effect of prehadronic interactions we show the modifications of the 
multiplicity ratio due to the conventional hadronic FSI after the formation time tj and 
explore the sensitivity of the results on the size of tj. We, therefore, neglect any interactions 
during tj and for simplicity assume that tj is a constant in the rest frame of each hadron 
independent of the particle species. Due to time dilatation the formation time tj in the 
laboratory frame is then proportional to the particle's energy (cf. Eq. (j!TT|) ). Fig. [3 shows 
the result for the ^^Kr target using formation times in the range tj = 0-1.5 fm/c. For Tf=0 
(dashed line) we get a much too strong attenuation both in the z^ and z/ spectrum. In 
this case all reaction products start to interact immediately after the 7*A^ interaction and 
there is no effect of time dilatation {tj = tj = 0). As a consequence this limit leads to an 
almost flat z/ dependence of R'l^. We note that without the cut on the hadron energy E^ 
and without the limitation by the HERMES detector acceptance one would find a strong 
increase of the hadron multiplicity at low Zh due to particle creation in the FSI (cf. Fig. Iiup. 
However, after applying all cuts one is left with a 10% effect only, as can be seen from the 
Zh spectrum in Fig. [7| 

A slight increase of the formation time tj from zero to 0.1 fm/c (dotted line) already leads 
to a dramatic change in R'lj. By using this unphysically small formation time one obtains 
a good description of the z/ dependence, but fails to reproduce the high energy part of the 
spectrum in Zh- The reason is that many of the high-energy particles, which are directly 
created in the primary 7*A^ interaction, now escape the nucleus due to time dilatation. 
Especially the formation times tf of the light pions start to exceed the dimension of the ^^Kr 
nucleus for energies larger than E-j^ ^ 13 GeV. One also sees that the multiplicity ratio at 
low Zh drops below one because both the absolute number of FSI and the energy available 
for particle production in the FSI is reduced. The latter decreases the probability that the 
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secondaries survive the experimental cuts. 

At larger values of Tf one still has some attenuation due to the FSI of more massive 
hadrons. As seen in Fig. El a large fraction of the finally detected pions stem from the 
decay of neutral (and charged) p mesons created in the string fragmentation. Due to their 
relatively large mass the latter have a smaller formation time in the lab frame and are subject 
to hadronic FSI, especially if they carry only a small fraction Zh of the photon energy. For 
example, the formation length of a p meson produced by a 7 GeV photon is always smaller 
than about 7 ■ r^ • c ~ 9 • rj ■ c. For tj = 1.5 fm/c (short dashed line) only the formation 
lengths of the massive vector mesons and nucleons with Zh < 0.6-0.7 are short enough to 
give rise to attenuation. We mention that the small deviation from unity of R^.^ at ^/^ ~ 1 is 
due to the Fermi motion of the bound nucleons that affects the maximum energy available 
for hadron production in the initial 7*A^ interaction. 

B. Constituent quark model 

From Fig. [7|one extracts that for reasonable formation times Tf > 0.5 fm/c the (pre-) 
hadronic interactions have to set in quite early after the 7*A^ interaction, especially for the 
hadrons at large z^. The Lund model predicts a vanishing prehadron production time as 
z/i — i> 1 as indicated by Eq. (I2 8|l. However, as we have shown in Ref. [uj the Lund production 
time of Bialas and Gyulassy |2S| does not reproduce the charged hadron multiplicity ratios 
if one also accounts for the production of secondary particles in the FSI. As in our previous 



approach 11[ we, therefore, set the production time of all prehadrons to zero, but now 
rescale the prehadronic cross sections according to the constituent quark concept ((221) • 

Fig. ISl shows the results of our simulation for ^^N and '^^Kr using formation times Tf = 0- 
1.5 fm/c. The dashed lines (rj = 0) coincide with the ones in Fig. 13 since they only involve 
hadronic FSI. By comparing the u spectra of Figs. |H1 and 13 for finite formation times one 
observes that R\,j{v) is reduced by the prehadronic FSI, which also improve the agreement 
of R'lj{zh) with the experimental ^^Kr data at large z^. However, the attenuation of the 
Zh spectra is too strong for the ^"^N target which might already indicate a deficiency of our 
simple prehadron concept. The experimental ^^Kr data favor formation times Tf >0.3 fm/c 
with only little sensitivity for larger values because of the finite size of the ^^Kr nucleus. 
These times are in line with our simple estimate via the hadronic radius and with the values 
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Tf = 0.4-0.8 fm/c 42| extracted from antiproton attenuation in pA reactions at AGS energies 



of 12.3 GeV and 17.5 GeV on various nuclear targets |4^. A cleaner discrimination between 
Tf = 0.3 fm/c and larger formation times will be possible in an experimental investigation 
at Jefferson Lab J44l] with heavier targets and lower photon energies. 

In Fig. El we show the transverse momentum dependence of the multiplicity ratio (|33|). 
where the transverse component pt of the hadron is defined with respect to the momentum 
direction of the virtual photon. In the simulation we use the constituent quark concept 
()32|) for the prehadronic cross sections and the formation time Tf = 0.5 fm/c. We expect 
that the px distribution of the observed hadrons is broadened for complex nuclei compared 
to deuterium due to multiple scattering of the (pre-)hadrons. Up to p^ ~ 1 GeV^ the p"^ 
dependence of the multiplicity ratio is well reproduced for both the nitrogen and the krypton 
target. However, the data of Ref. [3] show a strong increase of -R^j for Pj^ ^ 1 GeV^, which 
is not reproduced by our (pre-)hadronic FSI even if one assumes that all elastic scattering 
events are isotropic in the center-of-mass system (dashed line). This can be considered a 
signal for a partonic origin of the enhancement of high-pT hadrons in eA collisions either 
via a change of the parton distributions inside the nuclear medium and/or the Cronin effect 



1451 



46|,|47 



The Cronin effect was first observed in 1975 



via an enhancement of high-pT hadrons 



in pA collisions and has become especially important recently in connection with data from 



% Isf. Similar to pA colhsions |4fil |. 



a 



high-pT hadron production in heavy- ion collisions 

high-energy parton produced in a direct 7*A^ interaction may be subject to soft coherent 

and incoherent multiple rescatterings in the nuclear medium. While the incoherent rescatter- 

n 

ings can be interpreted intuitively as a random walk in transverse momentum space [48|, the 
coherent gluon radiation from different nucleons is subject to Landau-Pomeranchuk-Migdal 
interference effects |8j. The authors of Ref. [l3| have calculated the transverse momentum 
broadening that is caused by the multiple scattering of the struck parton before the pre- 
hadron production time. Their calculation reproduces the strong increase of R\j{pt) for 
Pt > 1 GeV/c and thus supports the interpretation in terms of a Cronin effect. 

In our simulation we do not find any dependence of the charged hadron multiplicity ratio 
(I33|) on the photon virtuality Q^ . This is no surprise since in the ratio R\j the Q^ dependence 
of the primary electroproduction cross section cancels out and our prehadronic cross sections 
(I32|l do not depend on the virtuality of the photon. Therefore, an experimentally observed 
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enhancement of R\j with Q^ could be interpreted as a signature for color transparency 



49, 



50|. We will show in Sec. IVIII that indeed the simple constituent quark ansatz for the 
prehadronic cross sections (|32j) overestimates the attenuation in the kinematic regime of the 
EMC experiment, i.e. at larger values of v and Q^. 

C. Acceptance cuts 

We now discuss how the geometrical acceptance of the HERMES detector and the kine- 
matic cuts affect the multiplicity ratio (j33|) . In Fig. ^J we compare the results of our 
simulation using the constituent quark concept (jH^ and a formation time r/ = 0.5 fni/c for 
the HERMES acceptance (solid line) in comparison to a 47r-detector (dashed line). In both 
calculations we still account for all kinematic cuts in the HERMES experiment. As can be 
seen from the z^ spectrum, a detector with full angular coverage (dashed line) will detect 
many more of the low-energy particles - produced in the FSI - which simply do not end 
up in the HERMES detector. As a result, the v spectrum for a 47r-detector is almost flat 
since an increase of the formation time with v due to time dilatation not only reduces the 
attenuation but also the particle production in the FSI. According to our simulations the 
slope in the v spectrum experimentally observed at HERMES partly arises because at lower 
photon energies particle production in the FSI is less forward peaked and, therefore, less 



particles are seen by the H 
other approaches 



IRMES detector. Note, that this problem is usually neglected in 



. m lid, Il3{ that intend to describe the observed multiplicity ratios. The 
dotted line in Fig. ^J shows the result of a simulation where in addition to the geometrical 
acceptance the Eh> lA GeV cut has been neglected. As can be seen from the low z^ part of 
the multiplicity ratio, one now detects even more \ow-Zh hadrons. Without the cut Zh > 0.2 
the multiplicity ratio as a function of the photon energy rises to R\,^{v) ~ 1.5. 

Summarizing this paragraph we stress that one has to take the geometrical acceptance 
of the HERMES detector into account if one wants to draw conclusions from a comparison 
of experimental with theoretical results for the Xow-Zh part and the z/ dependence of the 
multiphcity ratio. 
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D. Double hadron attenuation 

Before we turn to the individual attenuation of the various identified hadrons we compare 
ou. ..uMio„ With the .ecently .eas.ed doub.e-had..„n atte„uat,o. at HERMES g In 

each event only the two (charged or neutral) hadrons with the highest energies are consid- 
ered. In the following we denote the hadron with the highest Zh as the leading hadron and 
the other one as the subleading hadron. The experimental observable is the double-hadron 
attenuation ratio 

R2{Z2) = -j^^trr- (34) 



N2iz2) 


Ni A 


N2{Z2) 


Ni D 



Here A^2(-22) denotes the number of events where the leading and subleading hadron carry 
the energy fraction zi > 0.5 and Z2 < zi, respectively, and A^^i is the number of events where 
at least one of them has Zh > 0.5. The kinematic cuts are the same as for charged hadrons 
except for the Bjorken variable x, which now has the new boundary x > 0.01. 

Fig. ^2 shows the double-hadron multiplicity ratio ()34|) for ^^N, ^^Kr and ^^^Xe. To 
exclude contributions from p° decay into tt+tt^ the charge combinations 'H — ' and ' — h' 
have been excluded both in experiment and in the simulation. The solid line shows the 
result of a full coupled-channel calculation using the constituent quark concept (jH^ and the 
formation time Tf = 0.5 fm/c. The shape of the spectrum is similar to that of R\j{zh) of 
the charged hadron multiplicity ratio shown in Fig. |H| The reason is quite simple: For the 
interpretation we discard for a moment the constant factors A^^i in Eq. (jM|) and the factors A'^e 
in Eq. ()33|) which have no influence on the shape of the Zh dependence. The only difference 
between Nh{zh)\A/Nh{zh)\D and A^2(-Z2)U/A^2(-22)|d then is that one restricts the detected 
hadron to the subleading particle in the latter case. If the subleading particle (with energy 
fraction Z2) of the initial 7*A^ reaction interacts with the nuclear environment, it will produce 
a bunch of low-energy particles. The new subleading hadron in the event then has a energy 
fraction z'2 < Z2- As for the usual charged hadron multiplicity spectrum the coupled-channel 
FSI shuffle strength from the high z^ part to the low z^ part of the spectrum. This is not the 
case for purely absorptive FSI (dashed line in Fig. ITT|) . As one can see, our coupled-channel 
calculations (solid line) - using the constituent quark concept ()32|) and the formation time 
Tf = 0.5 fm/c - are again in quantitative agreement with the experimental data apart from 
the last data point in the ^^Kr data, which indicates a multiplicity ratio R2{z2 = 0.5) ~ 1. 
This behavior cannot be explained within our model. 
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As can be seen from Fig. [TTJ our calculations predict about the same double-hadron 
attenuation ratio for ^^^Xe and ^^Kr. The reason is that the attenuation of leading and 
double hadrons increases in the same way when going from the krypton to the xenon target. 
Hence, the double-hadron attenuation ratio (j34j) stays roughly the same. Note, that this 
does not necessarily imply the same hadron attenuation for *^Kr and ^^^Xe. As we have 
shown in Ref. [ll|, and as one can also see by comparing the multiplicity ratios for the two 
targets in Figs. ^1 and El the hadron attenuation in the ^^^Xe nucleus is on average 5% 
larger than for ^^Kr. 

E. Attenuation of identified hadrons 

We finally consider the attenuation of vr^, 7r°, i^^, p and p in DIS of 27.6 GeV positrons 
off ^°Ne, ^^Kr and ^^^Xe nuclei. For the krypton and xenon target the cuts are the same as for 
charged hadrons plus the momentum cuts necessary for particle identification at HERMES 
ja], i.e. 2.5 GeV/c < p^,/^ < 15 GeV/c and 4 GeV/c < Pp^p < 15 GeV/c. The following cuts 



are different for the neon target jsa]: 2 GeV< z/ < 24 GeV, 0.6 GeV/c < p^ < 15 GeV/c 
and 2 GeV/c < pk < 15 GeV/c. 

For the moment we stick again to the prehadron concept ()32j) and the formation time 
Tf = 0.5 fm/c in our simulation. Figs. IT^ and IT^ show the multiplicity ratios on '^'^Kr and 
^°Ne, respectively. The attenuation of pions is well described for the ^^Kr nucleus while it 
is slightly too strong for the light ^°Ne target. This is no surprise since for ^"^Kr our model 
already reproduced the multiplicity ratio of charged hadrons - dominated by pions - while 
the attenuation was too strong in the case of the lighter ^^N nucleus. 

The attenuation of K'^ mesons is well described for the heavier krypton target while it is 
in poor agreement with the ^°Ne data. In the latter case the calculated multiplicity ratio is 
larger than one at small Zh in contradiction to the experimental data. For both nuclei our 
simulation yields approximately the same attenuation for K~ and K'^ mesons at large Zh- 
The reason for this is directly related to the K^ and K~ spectra in Fig. which show that 
due to the quark content {us) the K^ contain less quarks from the beam or target than the 
K^. The few K^ that are not solely made of quarks and antiquarks created from the vacuum 
in the string fragmentation carry (anti-)quarks from the resolved photon component or the 
nucleon sea. According to the constituent quark concept (jH^ one has thus less prehadronic 
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interactions of K~ . This is compensated by the larger /^"-nucleon cross section. In total 
this leads to a similar attenuation of K^ and K~ . From Fig. El it can also be seen that 
a large part of the kaons at high Zh stem from decay into K'^K". The attenuation of 
K^ at high z^, therefore, strongly depends on the FSI of the (f) meson; this is in analogy to 
the attenuation of charged hadrons (pions) which are strongly affected by the FSI of the p 
meson. 

There is a further complication connected with the multiplicity ratio of kaons, which 
is neglected in a purely absorptive treatment of the FSI. The initial 7*A^ interaction pro- 
duces many more pions and p mesons than strange particles (cf. Fig. E)). These high-energy 
particles can produce secondary K^ and K^ in the nuclear FSI and thereby enhance the 
multiplicity ratio for K^ at low Zh- This is illustrated by the dotted lines in Fig. 1131 which 
shows the result of a purely absorptive treatment of the FSI, in which every particle that 
undergoes FSI is simply removed from the simulation. One clearly sees that kaon absorption 
in the FSI is compensated to a large extent by the production of kaons in the nuclear FSI 
of pions and p mesons. For Zh ^0.35 the K~^ production exceeds the absorption in the light 
^°Ne and leads to a multiplicity ratio larger than one. This is not the case for the *^Kr 
nucleus which is large enough to also absorb some of the secondary kaons. Of course this 
effect strongly depends on the strangeness production cross section used in the FSI. Unless 
one does not have all these coupled channels under control it is, therefore, hard to draw any 
conclusion from the kaon attenuation in DIS off nuclei. 

From Fig. El one also sees that there are only very few protons that contain no quarks 
from the beam or target because diquark-antidiquark creation is strongly suppressed in the 
string fragmentation due to the relatively large diquark masses. The latter also explains why 
most protons contain two of the original quarks, i.e. the diquark from the target nucleon. 
These remnants have a very large prehadronic cross section according to the constituent 
quark concept (J32|) . As a result the proton multiplicity ratio comes out too small for both 
the ^"^Kr and the ^°Ne nucleus. As discussed in Sec. IV A| the Lund model points to a larger 
formation time for protons than for the lighter pions, p mesons and kaons (cf. Eq. (jSHI))- 
However, as we have shown in Fig. jHl a further increase of the formation time to r > 0.5 
fm/c does not change our result at large Zh as long as the production time is zero. Since 
proton attenuation at large Zh is solely due to prehadronic interactions this either points 
towards a smaller prebaryonic cross section or a finite production time. 
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A further problem connected with the proton spectra is the strong increase of the mul- 
tiphcity ratio for Zh < 0.4 which is seen in the experimental '^^Kr data. At small proton 
energies i?^ becomes larger than one which might be understood in our model by a slowing 
down of high-energy protons in the FSI. Alternatively, protons might be knocked out of 
the nucleus in the FSI of a high-energy meson produced in the primary 7*A^ interaction. 
We do indeed see these effects in our simulation, however, in both cases the experimental 
momentum cut pp > A GeV/c removes most of these protons from the acceptance. Thus, the 
protons in our transport simulation loose too much energy per collision in the (pre-)hadronic 
FSI scenario. 

As we have pointed out in the discussion of Fig. |21the use of FRITIOF in the simulation 
of 7*A^ events has a small effect on the proton spectra. In Fig. ^J we, therefore, also show 
the proton attenuation in a simulation where all '-f*N interactions above Wpy = 3 GeV are 
simulated by PYTHIA. One observes that besides a slight improvement at large z^ and low 
photon energies z/ this has no effect on our result. 

Within our simulation the attenuation of antiprotons also comes out slightly too large 
in the (pre-)hadronic FSI scenario. We find that the antiprotons with Zh ^ 0.5 are mainly 
beam or target remnants that contain an antiquark from the resolved photon or the nucleon 
sea, whereas most of the antiprotons with Zh < 0.5 are solely made of antiquarks that are 
produced in the fragmentation of the string excited in the 7*A^ interaction. According to the 
constituent quark concept (J32|l the attenuation of antiprotons with Zh < 0.5 is, therefore, 
only caused by hadronic FSI after Tf. 

In Fig. El we also show predictions for the attenuation of identified hadrons on a ^^^Xe 
target. In the simulation we have used the constituent quark concept ()32|) and a formation 
time Tf = 0.5 fm/c. By comparing the multiplicity ratios R'l^ of negatively and positively 
charged pions for ^°Ne, ^^Kr, and ^^^Xe nuclei, we find that the attenuation (1 — -R^) 
scales like A" with an exponent a = 0.22-0.29 at Zh = 0.95. If the attenuation was simply 
proportional to the distance that the particles propagate through the medium, one would 
naively expect a scaling exponent a = 1/3. The deviation is caused by the finite formation 
time and the nonuniform density distribution in nuclei. At lower values of Zh or in the 
integrated z/ spectra the scaling behavior is hidden by the coupled channel effects. We note 
that a scaling of the attenuation with the target mass ~ (A)^/^, as predicted by Ref. 
would imply an increase of about 34% when using ^^^Xe instead of ^^Kr. 
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F. Lund production and formation times 

We now test the result of a simulation where both the production times of the prehadrons 
and the formation time Tf of each individual particle are explicitly extracted from the cor- 
responding string fragmentation in JETSET (cf. Fig. HJ. As described in Sec. IV Al there are 
in principle three time scales involved in the Lund fragmentation process: i) The produc- 
tion proper time Tpi of the hadron's first constituent, which is obviously zero if the hadron 
contains a constituent from a string end, ii) the production proper time Tp2 when the second 
constituent is produced and a color neutral object is formed, and iii) the formation proper 
time Tf where the two world lines of the constituents cross for the first time. 

We point out that - due to technical reasons - all particles that emerge from the primary 
7*A^ interaction start to propagate from the interaction vertex, while the production time 
only affects the beginning of their interactions with the nuclear medium. One may therefore 
expect slight deviations from the real reaction geometry, where the excited string propagates 
over a small distance prior to fragmentation. 

We now proceed and introduce two effective cross sections: Ui which accounts for the 
'partonic' interaction between r^i and Tp2, and 0"2 which accounts for the 'prehadronic' 
interactions of the color neutral object between Tp2 and Tf. We first neglect the 'partonic' 
interactions (cti = 0) and set the 'prehadronic' cross section a2 equal to the full hadronic 
cross section cxh. In such a scenario we are no longer sensitive to the formation time Tf 
but only on Tp2. The solid lines in Fig. ITHlshow the result of such a simulation for a ^^Kr 
target. Since the resulting attenuation is much too weak we conclude that the strong FSI 
have to start earlier. Note that a reduction of o"2 will further enhance the discrepancy with 
experimental data. To achieve reasonable results for the multiplicity ratio one has to rescale 
Tp2 by about a factor of 0.2 (dotted line), which is quite dramatic. According to the Lund 
model (cf. Eq. (j^Hj) ) this has to be interpreted as an increase of the string tension k, by an 
unreasonably large factor of about five in the nuclear medium. Setting ai to zero implies to 
neglect any interaction of the nucleon debris with the nuclear medium between the moment 
of the 7*A^ interaction and Tp2. This might be a problem, since the hadronic string that is 
produced in the DIS may interact with a hadronic cross section right from the beginning 



The average size of the prehadron-production times of the gluon-bremsstrahlung 
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model [l3[ is about a factor of ten smaller than the times tp2 extracted from JETSET. 
In fact, their average size is rather compatible with tpi. If one assumes that strong FSI 
already set in right after Tpi and sets cxi = a2 = (Jh one gets the result indicated by the 
dashed curves in Fig. ^[ which are in satisfactory agreement with the experimental data. 
However, such a large interaction cross section o"i is definitely not of perturbative nature. 
Furthermore, this recipe again implies that all beam and target remnants can interact right 
after the photon-nucleon interaction since the production time Tpi of their first constituent 
is zero. In fact, this scenario is not much different from our constituent quark ansatz dis- 
cussed before. Again the interactions of the beam and target remnants may be considered as 
effectively accounting for the interaction of the strings right after their creation in the 7*A^ 
interaction. The solid squares in Fig. E] indicate the average starting times of the prehadronic 
and hadronic interactions according to the constituent quark model ()32j) with Tf = 0.5. The 
Zh dependence of these two times has to be compared with tpi and tf, respectively. In both 
cases the shape looks quite similar while the average times of our constituent quark concept 
are somewhat smaller. The latter is partly compensated by the reduced prehadronic cross 
section, cf. Eq. fl32j) . Due to time dilatation the production times tpi are in general already 
of the order of the nuclear radius. This explains why the beam and target remnants - for 
which tpi = - dominate the shape of the spectra. 

VII. HADRON ATTENUATION AT EMC ENERGIES 

In this section we test different space-time pictures of hadronization in comparison to the 
EMC data with 100 and 200 GeV muon beams |lj]. In the previous section we have seen 
that almost all of the HERMES data can be described with the simple prehadron concept 
of Sec. IV B| i.e. setting the production time Tp to zero for all hadrons and using the con- 
stituent quark concept ()32|1 for the prehadronic cross sections during the formation time Tf. 
Obviously, this picture can only represent a rough approximation to the real hadronization 
process. Neither is it very likely that the string fragments convert instantaneously into color 
neutral prehadrons nor do the cross sections instantaneously iuiap from the rescaled values 
(I32|l to the full hadronic size. 

The kinematic regime of the EMC experiment, which uses a 100 GeV (200 GeV) muon 
beam, is different from that of the HERMES experiment. Here, the kinematic cuts are 
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g2 > 2 GeV^, W >4: GeV, x > 0.02, 10 GeV< z/ < 85 GeV (30 GeV< z/ < 170 GeV) and 
Eh > 3 GeV. In addition we again account for the angular acceptance of detector, i.e. ±5° 
horizontally and ±8° vertically [53[. 

In the upper panel of Fig.^Jwe show the EMC result for the multiplicity ratio of charged 
hadrons on a ^^Cu target as a function of the hadron energy fraction z^. The solid line shows 
the result for the simple prehadron concept ()32|) with the constant formation time tj. The 
shaded area indicates the region between the results of calculations using 100 GeV and 200 
GeV muon beam energy. Obviously, we get a much too strong attenuation for z^ > 0.3. This 
either implies that the prehadronic interactions set in too early or that the cross sections 
(IH^ are too large. 

The dashed line shows the result for a simulation when assuming that during Tf the 
prehadron cross section increases quadratically in proper time r from zero to the asymptotic 
value ah- Such an ansatz can be motivated by color transparency, which states that the 
cross sections of the color neutral prehadron scales with its diameter squared. While giving 
the right attenuation for EMC energies, such an ansatz fails to explain the HERMES data 
as long as the initial cross section is exactly zero (see lower panel of Fig. fTH|) . For the 
calculation indicated by the dash-dotted line in the lower panel of Fig. ^J we assumed an 
initial prehadronic cross section O.Sah followed by a quadratic increase in proper time up to 
the full hadronic value at- The difference between the initial prehadronic cross sections at 
HERMES and EMC might be explained if one assumes that the cross section right after the 
7*A^ interaction is set by the resolution of the virtual photon, i.e. Q^. Indeed, the average 
values of Q^ in the EMC experiment are more than twice as large as at HERMES energies. 



However, the HERMES data J5j] do not indicate a strong enough Q^ dependence of the 
multiplicity ratio to explain such a dramatic difference between the two 'initial' values for 
the prehadron cross section. 

In Sec. IV II we have shown that, as long as one neglects a strong partonic energy loss [s], |9| 
right after the 7*A^ interaction, the (pre-)hadronic interactions have to set in very early 
(cf. Fig. Ej) to explain the various HERMES data. Since the Lorentz 7 factors involved in 
the EMC experiment are about five times larger than at HERMES energies, a finite (but 
small) production time will have a larger impact on the calculated multiplicity ratios. In 
Fig- El we, therefore, also show the results of a calculation with a constant production time 
Tp = 0.1 fm/c as well as the production time Tp = 0.2x^2 with Tp2 extracted directly from 
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JETSET. For simplicity we neglect all (partonic) interactions before the production time of 
the prehadrons and set the prehadronic cross to the full hadronic cross section ah- 

The result for Tp = 0.1 fm/c in the kinematic region of the HERMES experiment is 
shown in Fig. [7| and yields a too weak attenuation at large z^. An additional formation time 
with reduced cross sections would further enhance this discrepancy. On the other hand, 
the production time Tp = 0.1 fm/c is still too small to give the right attenuation at EMC 
energies as can be seen from the dashed line in Fig. El 

In Sec. IVI Fl we have found that using r^i as the prehadron production time yields a sat- 
isfactory description of the HERMES data. However, using r^i as the prehadron production 
time leads to the same problem observed for our constituent quark concept in Fig. E[ i.e. a 
too strong attenuation of the high-Zh hadrons. The reason is again the strong prehadronic 
interactions of the beam and target remnants right after the 7*A^ interaction. 

When setting the prehadron production time to 0.2x^2 (solid line in Fig. ITTj). we observe 
again a slightly too strong attenuation. However, this concept is in better agreement with 
the experimental data at EMC energies than our previous approaches. 

In summary, it does not seem to be possible to simultaneously describe both the HER- 
MES and the EMC data with (pre-) hadronic FSI only, as long as one does not account for 
additional effects like color transparency. 

VIII. SUMMARY AND OUTLOOK 

In this work we have presented a detailed investigation of hadron attenuation in deep- 
inelastic lepton scattering at HERMES and EMC energies. In the primary production 
process we do not only consider the direct interaction of the virtual photon 7* with a 
quark inside the nucleon (nucleus), but also diffractive and hard scatterings of the photon's 
hadronic components. In the latter case we account for coherence length effects that lead 
to nuclear shadowing in DIS of nuclei (cf. Section ITTT]) . Furthermore, our model description 
incorporates other nuclear effects like Fermi motion of the bound nucleons, Pauli blocking 
and nuclear binding. We determine the complete final state of each 7*A^ interaction using 
the event generator PYTHIA. 

For the space-time picture of the hadronization process we basically have considered two 
different scenarios. One is motivated by the Lund model and the other one is that generally 
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used in standard transport models. In the Lund model the production proper times of 
the hadron constituents and the formation proper time of the hadron are determined by the 
underlying fragmentation process and depend on energy and momentum of the hadron. Since 
PYTHIA is based on the Lund fragmentation scheme we can directly extract the production 
and formation times for each hadron from the corresponding fragmentation process. We 
have shown explicitly that the production and formation times in the Lund model show a 
non-trivial dependence on the mass of the hadron. 

On the other hand, the concept generally used in transport approaches is that the produc- 
tion time of the prehadrons is set to zero and the interaction probability is reduced during 
the formation time Tf. For simplicity the formation time is assumed to be a constant Tf in 
the rest frame of each hadron in order to work with a single parameter, only. 

The prehadronic final state interactions (FSI) between the production and the formation 
time and the hadronic FSI after Tf are described within a transport model which allows for 
a realistic coupled channel treatment beyond simple absorption mechanisms. We explicitly 
account for particle creation in the interactions of the primary reaction products emerging 
from the initial 7*A^ interaction. These secondary particles are found to strongly influence 
the low-energy part of the experimentally observed multiplicity ratios ()33|) . 

Furthermore, we have studied how the kinematic cuts and the finite detector acceptance 
influence the experimental observables. We find strong effects that have to be taken into 
account in any robust interpretation of the data. 

We have investigated the attenuation of charged and neutral pions, kaons, protons and 
anti-protons as a function of the fractional energy Zh, the photon energy u and the transverse 
momentum p^ as well as the double-hadron attenuation. While in the kinematic region of 
the HERMES experiment most phenomena can be attributed to (pre-) hadronic FSI, we find 
limitations of our 'standard' model for EMC energies and for large pt- The latter supports 
a partonic origin for the Cronin effect in electron-nucleus interactions. 

We have seen in Sec. lVIFl that a scenario in which any interactions before the production 
of the color neutral prehadrons at proper time Tp2 are neglected does not yield enough hadron 
attenuation. Instead, we had to assume strong interactions right after the production times 
Tpi of the hadron's first constituent. However, this does not necessarily imply that there is 
a freely interacting quark that sees strong FSI. As pointed out in Sec. IV Al t he longitudinal 
dimension of the string and its (color neutral) fragments is not expected to be larger than 
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1 fm J36|, |37|. Therefore, one has in principle to deal with the propagation of color neutral 
strings that subsequently fragment into prehadrons and color neutral remainder strings. 
A collision of a (remainder) string - before its fragmentation - will most likely lead to a 
different final state than an undisturbed decay. This is technically more involved but will 
be incorporated in our future work. Apart from extensions on the theoretical side new 
experimental data also at lower energies should help to clarify the problem of hadronization. 
The experiments that are currently performed at considerably smaller energies at Jefferson 
Lab - and which are planned after the upgrade to 12 GeV beam energy - will be more 
sensitive to the string fragmentation and hadron formation times since time dilatation effects 
are less pronounced. 

Furthermore, in a direct photon-nucleon interaction at very high Q^ and energy v the 
photon is expected to knock out a highly virtual point-like quark that immediately may 
radiate gluons. The gluons then can split into quark-antiquark pairs and finally the various 
colored quarks and gluons combine to form the hadronic strings. While these processes are 
in principle taken into account in our present simulations via PYTHIA, we have neglected 
their space-time evolution so far and assumed that they take place instantaneously at the 
interaction point. This simplification may have considerable consequences at EMC energies 
and possibly explain our difficulties in describing the data. An explicit consideration of these 
partonic evolution effects is also planned for the future. 

Acknowledgments 

This work has been supported by BMBF. The authors acknowledge valuable discussions 
with A. Accardi, P. Di Nezza, G. Elbakyan, A. Hillenbrand, B. Z. Kopeliovich, P. Liebing, 
V. Muccifora, and X. N. Wang. 

APPENDIX A: HIGH-ENERGY HADRON-NUCLEON INTERACTIONS 

1. Total cross sections 

As pointed out in Sec. II VI we model hadron-hadron collisions above the resonance region 
[W > 2.3 GeV) by FRITIOF [23|. In contrast to PYTHIA |2lj the FRITIOF event generator 
does not provide absolute cross sections but only determines the final state of a scattering 

31 



event. Therefore, we have to exphcitly parameterize the total cross sections at large energies. 
For all baryon-baryon collisions - that solely involve non-strange baryons in the entrance 
channel - we use the cross section parameterization J55| 

a= [A + Bp'^ + C\n'^p + D\np]mh (Al) 

with the parameters for proton-proton collisions listed in Table Q] In Eq. ()A1|) p denotes the 



laboratory momentum in GeV. As shown in Ref. J55| this parameterization yields a good 
description of the cross sections down to the resonance region. In the current version of our 
transport simulation we neglect collisions of strange and charmed baryons. 

For the meson-nucleon interactions, that can be addressed experimentally like vr^p, K'^p 
and K^n, we use Eq. (jAlll with the parameters listed in Table HI The cross sections for the 
other isospin channels are given by isospin symmetry 

1, , 

'-'TT+n '-'■K~p 1 

'-'■K~n '^TT+p ) 

(TKOp = <^K+n , 

aaop = (TK-n, etc. (A2) 

The cross sections for the K* and K* mesons are assumed to be the same as for K and K. 
The total cross sections for the vector mesons are given by the PYTHIA parameterization 
©. 

For all other high-energy meson-baryon collisions we employ an ansatz similar to that of 
the UrQMD model |25]: 

^™1(^) = ^TnIs) (1 - OAx^s' - 0-5x^) (1 - OAx's - 0.5^), (A3) 

where x'g'^ and x^ q denote the strangeness and charm content of the meson m and baryon 

b, respectively: 

^m_\S{m)\ ^m_ \C{m)\ 

■^ ~ 2 ' "^ ~ 2 ' 

Here S and C denote the strangeness and charm quantum numbers of the meson and baryon, 
respectively. 
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2. Elastic scattering 

There are two problems with the original FRITIOF model 23|. The first one is that it does 
not generate enough elastic scattering events as can be seen by the solid triangles in Fig. ^1 
We cure this deficiency by also parameterizing the elastic cross section and simulating elastic 



scattering externally; FRITIOF is then only called for inelastic collisions (cf. Ref. J25|)- ^'^^ 
all elastic baryon-baryon collisions, that solely involve non-strange baryons in the entrance 
channel, we use the same parameterization as for elastic pp scattering (see Table |l}. 

The cross sections for elastic n^p and K^p scattering are again taken from experiment 
(see Table EJ) and those for elastic V^A^ scattering are given by VMD, 

el f9v\^ 



For the remaining particles we use the relation 

c^el = co■f/t^ (A5) 

to determine the size of the elastic cross sections from the corresponding total cross sec- 
tions. The value c = 0.039 mb^^/^ is the same as in the UrQMD model [35 1. The angular 
distribution for elastic scattering 1 + 2— i>l + 2is taken from the PYTHIA model. 



d& 



el 



dt 



~e^i2*, (A6) 



with the slope parameter 



Bf^ = 2bi + 262 + 4s' - 4.2 



with s given in units of GeV^ and Bf2 in GeV ^. The constants are 65 = 2.3 for baryons 
and 6m = 1-4 for mesons. 

3. Quark- ant iquark annihilation 

A second problem with FRITIOF is that it does not account for quark- ant iquark annihi- 
lation in meson-baryon scattering in correspondence to Reggeon exchange in the t channel. 
This process, however, gives a finite contribution to the total cross section for low ^/s. We, 
therefore, keep the option to simulate its contribution to the inelastic collisions indepen- 
dently of FRITIOF using the method developed in Ref. J26|, where the antiquark from the 
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meson and a constituent of the baryon with the same flavor may annihilate. We assume 
that the momenta of the two annihilating quarks are very small so that we do not need to 
treat the final gluon explicitly. The final state of such an annihilation process is modeled by 
an excited string with the invariant mass of the colliding system. The decay of the string 
into hadrons is taken into account using the Lund fragmentation routine JETSET 7.3 as 
also used by FRITIOF. In contrast to the two excited strings of a FRITIOF event, the 
string - emerging after qq annihilation - has a larger invariant mass and, therefore, has 
more energy available for ss creation. Hence, the annihilation process has a strong effect 
on strangeness production in meson-baryon scattering and the energy dependence of its 
cross section. Above ^ = 2.2 GeV the annihilation cross section relative to strangeness 
production in pion-nucleon scattering can be fitted by |2fil |: 



"''' -max 



^inel 



1.2-0.2^,0 



Fig. ^1 shows the resulting cross section for strangeness production in n^p with (solid line) 
and without (dashed line) the quark-antiquark annihilation contribution. One clearly sees 
that the incorporation of the annihilation part leads to a much better agreement with data. 
Since the cross sections for baryon-baryon collisions below ^/s = 2.6 GeV and meson- 
baryon collisions below y/s = 2.2 GeV are given by the resonance model of Ref. |l6|, we 
continuously connect the high-energy parameterization to the resonance part. As an exam- 
ple. Fig. ^1 shows our parameterization of the total and elastic n^p and K^p cross section 
in comparison with the experimental data. 

4. Antibaryons 

Finally, there is the possibility of elastic and inelastic baryon-antibaryon (bb) interactions 
in the nuclear FSl. The total and elastic cross sections for pp scattering are again taken from 
experiment, i.e. Eq. (jAlj) and Table|l] For the rest of the non-strange antibaryons we use the 
same cross sections as for antiprotons. Elastic bb scattering is simulated in the same way as 
in all other elastic channels. The inelastic fraction of the total cross section is experimentally 
known to be dominated by annihilation. We, therefore, reduce the inelastic cross section of 
(anti-)baryons that involve {s)s quarks according to the simple valence quark picture: 
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where 5* denotes the strangeness of the (anti-)hyperons. The annihilation in colhsions in- 
volving charmed (anti-)barvons is neglected here. The annihilation process is modeled in the 
same way as in Ref . J24|] : After annihilation of a quark and an antiquark with the same flavor 
the remaining (anti-) quarks form two orthogonal qq jets which equally share the invariant 
mass of the colliding system. As in the case of qq annihilation in meson-baryon scattering 
the strings are fragmented using JETSET 7.3. 

Throughout this work we neglect meson-meson, meson-antibaryon as well as antibaryon- 
antibaryon interactions. They are very unlikely in lepton-induced reactions since they require 
interactions between the few reaction products among each other. 

APPENDIX B: THE BUU TRANSPORT MODEL 

1. Ingredients 

The scalar potential Us in the set of transport equations (fTUj) can be related to an effective 
(non-relativistic) potential U which accounts for the many-body interactions of the baryons 
among each other. The general expression for the relativistic energy of a particle under the 
influence of a scalar potential S and a vector potential V = (Vq, V) is 



H=^{fi + Sy + (p- y)2 + Vo. (Bl) 

In the local rest frame (LRF), i.e. where the baryon current locally vanishes, the spatial 
components V vanish. We arbitrarily set 5* = and interpret the effective potential U as 
the zeroth component Vq of the vector potential: 



^LRF = ^J^i'' + pIrf + U{r, Plrf)- (B2) 

We can then define the scalar potential Us of Eq. (fTTj) in any frame as 



^5:=V^lrf-Plrf-/^- (B3) 

Note that for photon- and electron-induced reactions the local rest frame coincides with the 
frame where the target nucleus is at rest, i.e. the laboratory frame in case of fixed target 
experiments. 

For nucleons the effective potential U is parameterized according to Ref. [59] as a sum of 
a Skyrme part, which only depends on the baryon density p, and a momentum dependent 
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part: 

Cy.fl ^A^^B (m\ ^, f^ f^'/\ (B4) 

Po V Po 7 Po J (2t)''i I /p-p') 

where po = 0.168 fni~^ denotes the saturation density of nuclear matter. In the reactions 
considered in this work the nucleus remains close to its ground state and the phase-space 
density / in ()B4|1 can be approximated by the phase space density of (uncorrelated) cold 
nuclear matter 

/(r,p-)~e(p^(r)-|p-l) (B5) 

with the local Fermi momentum 

r 2 \ 1/3 

67r^ ^ ' 



Here g = 4: again denotes the factor of degeneracy. For the density distribution of complex 
nuclei we use the Woods-Saxon parameterization: 

1 + exp [—^) 
with the parameters of Table IVl that have been extracted from a Hartree-Fock calculation 
60| for stable nuclei. For light nuclei like ^^N, however, we use a Gaussian shape: 

P.W = ^exp(--^) (B8) 

with QgA^^^ = A/frrms and rrms = 1.21y4^/^ fm. 

The use of (IB5|1 allows us to employ an analytic expression for the momentum dependent 
part of the potential ()B4|1 



2C r dy eipFJr) - \p'\) ^2C g . 3 / pM^) + A' - p' Up + pFif)Y + A^ ' 
PoV (27r)3 i+fp^y ~ Poi2nr'' { 2pA ^ [(^ _ p^(f))2 + A2_ 

(B9) 



p + PF{r)\ fp-PF{r) 

arctan : — arctan 



2p^_ 
+ A \ ^ A J V A 

The parameters of the mean-field potential ()B4|) are fitted to the saturation properties of 
nuclear matter and the momentum dependence of the nucleon optical potential as measured 
in pA collisions [IG]. Throughout this work we employ the parameters given in Table ITTl 

We use the same mean-field potential for all baryons except for the A resonance for which 
we assume 

Ua = lu. (BIO) 
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This choice is motivated by the phenomenological value of -30 MeV at density po [SI]- As 
mentioned before, we neglect any hadronic potential for mesons as well as any influence of 
the Coulomb potential in our present investigations. 

In the multi-GeV range of interest in this work the coupling of the BUU equations (jlUp 
via the mean field is rather low because the absolute value of Us is about 100 MeV at 
maximum. One is, therefore, left with the coupling through the collision term (fT^ . In 
case of high-energy lepton- induced reactions binary collisions ai,a2 — * bi,...bm play the 
dominant role. Since some of the particles in our simulation are unstable with respect to 
strong decays one has to keep in mind that also the decay of a particle into several other 
hadrons leads to changes in the phase-space densities. 

In the mesonic sector we account for all particles listed in Table UTTl Mesons that can only 
decay due to the weak interaction are considered to be stable within our model. The only 
exception is the r] whose decay is explicitly accounted for when evaluating pion production. 

Besides the nucleon (m^v = 938 MeV) and the Delta {niA = 1232 MeV, Ta = 118 MeV) 
we account for 29 other nucleon resonances 1^ in the baryonic sector whose properties 
are taken from an analysis of vrA^ scattering |62|. In our model the A always decays to 
Ntc whereas the other resonances can couple to the channels A^vr, Nrj, AK, Nuj, Avr, Np, 
Na, iV(1440)7r and Ap. In addition to the A (mA = 1116 MeV) and the S (ms = 1189 
MeV), which are stable with respect to the strong decay, we also include 19 further S = —1 
resonances that can decay into Avr, NK, Svr, S*7r, Xrj, NK* and A*7r. Furthermore, we 
include the strange and charmed baryons of Table IIVI in our model. Due to a lack of a 
complete analysis, the parameters for the strange and charmed baryons are taken from 
Ref. J30|. For each baryon we also account for the corresponding antiparticle. 

Since we explicitly consider the charge of the particles each isospin state of a particle 
leads to a separate BUU transport equation. The spin is only accounted for as a statistical 
weight in the degeneracy factor g. 
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2. Numerical realization 

The set of coupled differential-integral equations (fTUj) is solved via a test-particle ansatz 
for the spectral phase-space densities (|12|1: 

1 (27t)^ ^ 
F{f,P,fi;t) = ^^^<5(f-rl(t))5(p-p,(t))(5(^-^,(t)) (Bll) 

N 9 tr 

where r^, pi and /ij denote the position, momentum and mass of the test particle i at time t 
and A^ is the number of test particles per physical particle. In this work we use the method 
of parallel ensembles, i.e. the test particles are divided into A^ different ensembles which do 
not influence each other. This is equivalent to simulating A^ independent nuclear reactions 
in parallel and averaging the observables at the end. For a test-particle number N —>■ oo 
the test particles will give the time evolution of the spectral phase-space densities. 

When initializing a nuclear reaction the test particles, that correspond to nucleons of the 
nucleus, are distributed in position space following a Woods-Saxon distribution ()B7|) or a 
Gauss distribution ()B8|) for ^^N, respectively. We here assume that the form of the density 
distribution is the same for protons and neutrons. For the initialization in momentum space 
we use the local Thomas- Fermi approximation ()B5jl . 

The calculation is performed on a discretized time grid with a default grid size At = 
0.5 fm/c. During each time step the test particles are assumed to move as non-interacting 
particles in the mean field Us- Substituting the test-particle ansatz ()B11|) into the BUU 
equation (fTUj) - with the collision term set to zero - yields the classical Hamilton equations 
of motion 



dfi 
dt 


= %.H 


dpi 

dt 


= -^nH 


d^i 
dt 


= 0, 



(B12) 

with H being a functional of the phase-space density /. 

Between the time steps the particles may collide. We do not assume any medium mod- 
ification of the matrix elements Ai that enter the collision term (fT^ . If one accounts for 
the energy shift caused by the scalar potential Us, the transition rates can be directly taken 
from the corresponding vacuum cross sections. Note, however, that in the resonance region 

38 



cross sections might be modified due to in-medium changes of the resonance properties as 
discussed in Refs. lid hn . 



Concerning the colhsion criteria we follow the method by Kodama et al. J63[: Two 
particles collide in a time step At if the impact parameter b, i.e. the minimum separation in 
their center-of-momentum system, is smaller than 



cr^°Hs) 



h<\l—^. (B13) 

Furthermore, it is checked if both particles reach this minimal distance during the time step 
At. In Eq. ()B13|1 a^"^ denotes the total cross section for the interaction of the two particles. 
For high-energy collisions these are the ones given in Appendix El Elastic interactions occur 
with a probability 

^el = 4^y (B14) 

the scattering angle is determined according to Eq. ()A6|) . In case of a high energy in- 
elastic collision the reaction products are determined by FRITIOF (or JETSET for baryon- 
antibaryon annihilation). In the resonance region, i.e. below ^/s = 2.2 GeV for meson-baryon 
and ^/s = 2.6 for baryon-baryon scattering, the total cross section is an incoherent sum of 
the cross sections for the reactions 

mB *-^ R, ttN ^ ttN, ttN -> vrvrA^, vrA^ ^ r/A, vrA^ ^ uN, ttN -^ ttujN, 

ttN ^ (pN, uN ^ uN, uN -^ tcttN, (pN ^ (pN, (pN -^ ttttN, tiB ^ KY, 

nB^KKN, KN^KN, KN ^ ttY, KN^nY*, KN ^ KN, KN ^ KnN 

in case of meson-baryon collisions and 

NN *^ NN, NN^NR, NN ^ AA, NN ^ NNn, NN ^ NNu, 
NR ^ NR', BB -^ NYK, BB -^ NY*K, BB -^ NNKK 

for baryon-baryon collisions. Here m stands for a meson, B=N, A; nucleon resonances 
are denoted by R and R', hyperon resonances by Y* and y=S,A. The reactions involving 
antibaryons are obtained by charge conjugation. The reaction channel ab ^ f in the collision 
of two particles a and b is chosen by Monte Carlo with a probability determined from its 
contribution to the summed total cross section (see Ref. |la| for details) 



ab^f 



CT, 



tot / 
ab ^ 



(B15) 
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An important feature of our model is the decay of unstable particles with mass /i and 
energy E during a time step At. The corresponding decay probability is given by 

r(^) 



Pdec = 1 - exp ( ^At j , (B16) 

where 7 = E/fi is the Lorentz factor while F denotes the width of the particle in its rest 
frame. The final state of the decay is again determined by Monte Carlo assuming the decay 
to be isotropic in the rest frame of the particle since we neglect the spin degree of freedom. 
Due to the low densities of other baryons Pauli blocking is only accounted for in collisions 
and decays that involve nucleons in the final state. For lepton-induced reactions, where the 
nucleus approximately stays in its ground state, one can approximate the probability that 
an event with an outgoing nucleon of momentum p is Pauli blocked via Eq. ()B5|) as 



Ppauli = e(p^(f) - \p\). (B17) 
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A 


B 


n 


C 


D 


4f 


48.0 


0. 


- 


0.522 


-4.51 


< 


11.9 


26.9 


-1.12 


0.169 


-1.85 


-%' 


38.4 


77.6 


-0.64 


0.26 


-1.2 


< 


10.2 


52.7 


-1.16 


0.125 


-1.28 


<% 


16.4 


19.3 


-0.42 


0.19 


0. 


<^p 


0. 


11.4 


-0.4 


0.079 


0. 


^tot 

■K p 


33.0 


14.0 


-1.36 


0.456 


-4.03 


■K-p 


1.76 


14.0 


-1.36 


0.456 


-4.03 


<+p 


18.1 


0. 


- 


0.26 


-1. 


4+p 


5. 


8.1 


-1.8 


0.16 


-1.3 


-.tot 


18.7 


0. 


- 


0.21 


-1.3 


tot 
^K-p 


32.1 


0. 


- 


0.66 


-5.6 


^K-p 


7.3 


0. 


- 


0.29 


-2.4 


tot 
^K-n 


25.2 


0. 


- 


0.38 


-2.9 



TAB 

Ref. 



E I: Parameters for the cross section parameterization (|A1|) . The parameters are taken from 

55|. 



A [MeV] 


B [MeV] 


C [MeV] 


T 


A [fm-i] 


-29.3 


57.2 


-63.5 


1.76 


2.13 



TABLE II: Parameters of the nucleon potential (|B4|) used in this work. 
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m 

[MeV] 


r 

[MeV] 


J 


/ 


s 


c 


decay channel 


vr 


138 








1 










V 


547 


1.2- IQ-^ 














77 (40%), TT+TT-TT^ (28%), 37r0 (32%) 


P 


770 


151 


1 


1 








vrvr 


a 


800 


800 














vrvr 


UJ 


782 


8.4 


1 











TTTT (2%), 7r07 (9%), vr+vr-^o (89%) 


i 


958 


0.2 














p% (31%), TTTTT] (69%) 


^ 


1020 


4.4 


1 











pvr (13%), KK (84%), tt+tt-tt^ (3%) 


K 


496 








1/2 


1 







K 


496 








1/2 


-1 







K* 


892 


50 


1 


1/2 


1 





Ktt 


K* 


892 


50 


1 


1/2 


-1 





Kvr 


Vc 


2980 



















J/^ 


3097 





1 













D 


1869 








1/2 









D 


1869 








1/2 





-1 




D* 


2007 


1 


1 


1/2 







Dtt 


D* 


2007 


1 


1 


1/2 





-1 


Dtt 


Ds 


1969 











1 






Ds 


1969 











-1 


-1 




Dl 


2112 


1.9 


1 





1 




Ds-i (94%), DsT^ (4%) 


D* 


2112 


1.9 


1 





-1 


-1 


Ds-i (94%), DsTT (4%) 



TABLE III: List of mesons included in the BUU model, their quantum numbers (spin J 



/, strangeness 5, charm C) and decay channels. The parameters are taken from Ref. 
for the a meson 



isospm 



30|] except 
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m 

[MeV] 


r 

[MeV] 


J 


I 


s 


C 


decay channel 


s 


1315 





1/2 


1/2 


-2 







"* 


1530 


9.5 


3/2 


1/2 


-2 





Svr 


n 


1672 





3/2 





-3 







Ac 


2285 





1/2 












Sc 


2455 





1/2 


1 









S* 


2530 


15 


3/2 


1 







Acvr 


^c 


2466 





1/2 


1/2 


-1 








2645 


4 


3/2 


1/2 


-1 




Hcvr 


f^c 


2704 





1/2 





-2 







TABLE IV: List of baryons with 5<— lorC>0 included in the BUU model, their quantum 
numbers (spin J, isospin /, strangeness S, charm C) and decay channels. The parameters are 
taken from Ref. 2i|- 





R [fm] 


a [fm] 


Po [fm-3] 


14n 


2.476 


0.479 


0.161 


20Ne 


2.851 


0.479 


0.161 


63 Cu 


4.409 


0.477 


0.157 


84Kr 


4.911 


0.476 


0.155 


i3iXe 


5.777 


0.476 


0.152 



TABLE V: Woods-Saxon parameters ()B7|1 for the nuclei investigated in this work. 
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FIG. 1: Parameterization of the cross sections for exclusive vector meson photoproduction ^p -^ Vp 



as a function of the invariant energy W in comparison to the data 



y. 



N 
T3 




0.2 0.4 0.6 



FIG. 2: Energy spectra of vr , K , p and p within the HERMES detector acceptance for a hydrogen 
target. The spectra are normahzed to the number of deep inelastically scattered leptons. The sohd 
hues represent the calculated results using our standard method of generating all resolved photon- 
nucleon events below Wpy = 4 GeV and all VMD events above Wpy with FRITIOF. Generating 
all events down to T^ = 3 GeV with PYTHIA leads to the spectra shown in terms of the dashed 
lines. The dotted lines show the result of a simulation without employing any kinematic or detector 
cuts on the hadrons. The data (full squares) are taken from Ref. jJi 
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FIG. 3: Simplified space-time picture of tlie fragmentation of a goQo string with invariant mass 
squared W^ = PqPq ■ The string decays into n hadrons with masses mi, . . . , rrin due to the creation 
of n — 1 quark-antiquark pairs from the vacuum. The production proper times Tp{j — 1) and Tp{j) 
correspond to the vertices where the Qj-iqj-i and the qjQj pairs are produced. The formation 
proper time Tf{j) of hadron j refers to the space-time point where the world lines of the quark 
Qj-i and the antiquark qj cross for the first time. Note that, the 'production times' of the qq and 
% - which in the end are contained in the first and the last rank hadrons - are always zero. 
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FIG. 4: Average (proper) times for hadron formation and the production of the first and second 
constituent as a function of z^ = Eh/v in the kinematic regime of the HERMES experiment. Ah 
proper times have been extracted directly from the JETSET routines which describe the string 
fragmentation in PYTHIA. The sohd squares indicate the average starting times of the prehadronic 
and hadronic interactions in the constituent quark concept (|32|) using a formation time rj = 0.5 
fm/c. 
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FIG. 5: Energy spectra of hadrons produced in electron-nucleon interactions at HERMES. The 
different lines indicate the hadrons that contain zero (sohd hne), one (dashed Une), two (dotted 
hne) or three quarks (dash-dotted hne) from the beam or target. For the proton the sohd hne 
nearly coincides with the z^ axis, which implies that most protons contain at least one leading 
quark. 
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FIG. 6: Model predictions for the average values of the kinematic variables in ch arg ed hadron 



production off ^^Kr in comparison with the experimental numbers at HERMES 4l|. For the 
calculation we used the formation time tj = 0.5 fm/c and the constituent quark concept (|32j) for 
the prehadronic cross sections. Left: {v) and (Q^) as a function of Zh compared to the experimental 
values for a ^"^Kr target (open symbols). Right: Same for (zh) and (Q^) as a function of v. 
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FIG. 7: Calculated multiplicity ratio of charged hadrons for the Kr target (at HERMES) assuming 
no prehadronic interactions during the formation time t/ = 0-1.5 fm/c. The data are taken from 
Ref. 
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FIG. 8: Calculated multiplicity ratio of charged hadrons for ^^N and ^^Kr nuclei (at HERMES) 
using the constituent quark model ()32() for the prehadronic cross section and different values of the 
formation time Tf = 0-1.5 fm/c. The data are taken from Ref. [4]. 
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FIG. 9: The multiplicity ratio of charged hadrons for ^"^N and Kr (at HERMES) as a function 
of the transverse momentum squared p^. In the simulation we use the constituent quark concept 
H32|) for the prehadronic cross sections and a formation time Tf = 0.5 fm/c. In the simulation 
indicated by the dashed line we additionally assume that_all elastic scatterings are isotropic in 
their center-of-mass frame. The data are taken from Ref. 
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FIG. 10: Influence of the HERMES detector geometry on the observed multiphcity ratios for a 
^^Kr target. In the simulation we use the constituent quark concept (|32() for the prehadronic cross 
sections and a formation time Tf = 0.5 fm/c. The sohd line represents the result of our simulation 
when accounting for the geometrical acceptance of the HERMES detector. In the simulation 
indicated by the dashed line no acceptance cuts have been employed. The dotted line represents 
the result of a simulation where in addition to the detector acceptance the Eh > 1.4 GeV cut has 



been neglected. The data are taken from Ref, 
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FIG. 11: Double-hadron attenuation ratio R2 for a ^^N, ^^Kr and ^^^Xe target as a function of the 
energy fraction Z2 of the subleading hadron. In the simulation (sohd hne) we use the constituent 
quark concept (|32j) for the prehadronic cross sections and a formation time tj = 0.5 fm/c. To 
exclude contributions from p^ decay into vr^vr^ the charge combinations 'H — ' and ' — +' have been 
excluded both in experiment and in the simulation. The dashed line shows a calculation with a 
purely absorptive treatment of the FSI. The preliminary HERMES data are taken from Ref. J5ll |. 
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FIG. 12: Multiplicity ratios of vr •'', K , p and p for a ^^Kr nucleus (at HERMES) as a function 
of the hadron energy fraction Zh = Eh/u and the photon energy v. The solid line represents the 
result of a simulation, where we use the constituent quark concept 1)32(1 for the prehadronic cross 
sections and a formation time Tf = 0.5 fm/c. The dotted line in the proton spectrum indicates 
the result of a simulation where all 7*A^ events are created by PYTHIA. The data are taken from 
Ref. y. 
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FIG. 13: Multiplicity ratios of vr^, K^, p and p for a ^'^Ne nucleus (at HERMES) as a function of 
the hadron energy fraction z^ and the photon energy z/. The solid line represents the result of a 
simulation where we use the constituent quark concept 1)32(1 for the prehadronic cross sections and 
a formation time tj = 0.5 fm/c. The dotted line represent the result of a simulation with a purely 
absorptive treatment of the FSI. The data are taken from Ref. |52l |. 
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FIG. 14: Predictions for the multiplicity ratios of vr , vr^, K , p and p for a ^^^Xe nucleus (at 
HERMES) as a function of the hadron energy fraction Zh and the photon energy z/. In the simulation 
we use the constituent quark concept 1)321) for the prehadronic cross sections and a formation time 
Tf = 0.5 fm/c. 
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FIG. 15: Multiplicity ratios of vr , vr*^, K , p and p for a ^^Kr nucleus (at HERMES) as a function 
of the hadron energy fraction Zh and the photon energy i'. In the simulation we use the proper 
times Tp2 (solid line), 0.2rp2 (dotted line) and Tp-^ (dashed line) from the JETSET routine as the 
prehadron production time. The prehadronic cross section is set to the full hadronic cross section 
and interactions before the production time are neglected. 



58 



1.2 
1.0 
0.8 
0.6 




a: 



1.0 



0.8 



0.6 



-1 ' r 



EMC 
■Standard BUU 

prehadron ^ f ' n 



H h 



H ' h 



^ 



HERMES 



■ a , , =min(T /T,,1)a 

prehadron ^ f ' h 

■CT ,, =min(0.3+TV,1W 

prehadron ^ f ' h 
J I I I I I I 



0.0 0.2 0.4 0.6 0.8 1.0 



FIG. 16: Upper panel: Multiplicity ratio of charged hadrons for a ^^Cu target as a function of 
the fractional energy Zh in the kinematic regime of the EMC experiment. The shaded areas are 
bounded by simulations using a 100 GeV (lower boundary) and 200 GeV (upper boundary) muon 
beam. The data are taken from Ref. jl^]. Lower panel: Multiplicity ratio of charged hadrons for 



a ° Kr target as a function of the fractional 



a 



energy z^ in the kinematic regime of the HERMES 



experiment. The data are taken from Ref. 121 • The solid line shows the result of our simulation 
using the constituent quark concept (|32() for the prehadronic cross sections and a formation time 
Tf = 0.5 fm/c. In the simulation represented by the dashed line we assumed a prehadronic cross 
section increasing quadratically in proper time during tj = 0.5 fm/c from zero up to the full 
hadronic size. The dash-dotted line represents a simulation where the prehadronic cross section 
increases quadratically in proper time from 0.3ah to the full hadronic cross section. 
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FIG. 17: Multiplicity ratio of charged hadrons for a ^^Cu target as a function of the fractional 
energy z^ in the kinematic regime of the EMC experiment. The shaded areas are bounded by 
simulations using a 100 GeV (lower boundary) and 200 GeV (upper boundary) muon beam. The 
data are taken from Ref . jlj | . The dashed line shows the result of our simulation using a constant 
prehadron production time Tp = 0.1 fm/c. In the calculation indicated by the solid line we used 
the Lund production time Tp = 0.2Tp2, where Tp2 has been directly extracted from JETSET. In 
both calculation we set the prehadronic cross section equal to the full hadronic cross section and 
neglect interactions before the production time. 
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FIG. 18: Cross sections for ir p and K p scattering. The solid (dashed) hnes represent the total 
(elastic) cross sections. The high-ener gy p art given by Eq. HA1|) is continuously connected to the 



cross sections of the resonance model 



id ] applied below ^/s = 2.2 GeV. The solid tr iang les show 



the FRITIOF result for elastic pion-proton scattering. The data are taken from Ref. 57 1. 
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FIG. 19: Cross sections for tt p — > strange particles in comparison with experimental data from 
Ref. Sg]. The solid (dashed) lines represent the simulation with (without) the possibility of quark- 



antiquark annihilation. The figure is taken from Ref. [2; 
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